o 



(N 






> 



Characterizations of distributions via order 
statistics with random exponential shifts 

M. Ahsanullah V.B. Nevzorov G.P. Yanev 



Abstract 

^\J ■ We characterize probability distributions via equalities in law between 

two order statistics shifted by independent exponential variables. An ex- 
plicit formula for the quintile function of the identified family of distri- 
butions is obtained. The results extend some known characterizations of 
exponential and logistic distributions. 
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1 Introduction 



Distributional (in law) relations between order statistics are useful and elegant 
tools for characterizing probability distributions. More specifically, there is a 
r>^ , large number of publications on characterizations based on recurrences involving 

Tlj" ' both order statistics and standard exponential variables. An excellent review 

r«^^ , and discussion of the available results on this subject can be found in [31. 

t^ ' AlZaid and Ahsanullah [5] showed that a non-negative random variable X 

with an absolutely continuous cumulative distribution function (cdf) has stan- 
dard exponential distribution iff for a fixed fc, such that 1 < fc < n — 1, 

KJl ; Xk,n + ;- = Xk+l,n, (1) 

C^ ' where Xi^„,X2,„ . . . ,X„^„ are the order statistics in a sample with parent X 

and ^ is standard exponential and independent from Xk,n- Please, see also [9] 
and [7j for an alternative proof and some extensions in the context of random 
contractions (cf. [6], p. 14). On the other hand, it was shown in Ahsanullah et 
al. [T], under some additional regularity assumptions, that a random variable 
X with absolutely continuous cdf is standard logistic if and only if for a fixed k 
such that 1 < k < n — 1, 

Xk,n H r = Xk+l,n T) (2) 

n — k k 

where ^i and ^2 are independent standard exponential variables, which are also 
independent from Xk,n and Xk+i,n, respectively. 



First, we will study some extensions of the distributional relations (H]) and ^ 
involving order statistics shifted by independent exponential random variables. 
What can be said about the distribution of the parent variable X if the following 
more general distributional equality holds 

Xk^n + a^i ^ Xr.n - b^2, a > 0, 6 > 0, I < k < r < n, (3) 

where ^i and ^2 are independent standard exponential variables, which are also 
independent from Xk^n and Xr^n, respectively? 

Our first result answers the above question in the case of adjacent order 
statistics, i.e., when r = fc + 1. Throughout this paper we will use the term 
Q{y) to refer to the pseudo-inverse (quintile) function of F{x), i.e., Q{y) = 
mi{x : F{x) > y} for y G (0, 1). 

Theorem 1 Let k be a fixed integer such that 1 < k < n — 1 and a > 

and b > be two fixed real numbers such that (a, 6) 7^ (0,0). Assume X is a 
random variable with continuous cdf F{x) and ^1 and £,2 are independent stan- 
dard exponential variables, which are also independent from Xk^n o,nd Xfc+i^„, 
respectively. Then 

Xk,n + a(l == ^fc+l,n - &6 (4) 

holds if and only if the quintile function Q{x) satisfies 

Q{x)^\ogLx''''{l-x)-'''^"-''A, 0<a;<l, (5) 

where c > is an arbitrary constant. 

Remarks. The following particular cases of (|H) and ([5]) might be of in- 
dependent interest, (i) Setting a = l/{n — k) and b — Q va ([S]), we obtain 
F{x) = 1 — e~'^~'^^ for x > c, i.e., an exponential cdf. (ii) If a = and 6 = 1 we 
have F{x) = exp{a; — c} for — cxo < x < c. (iii) If bk — a{n — fc), then ([S|) yields 
the logistic cdf (cf. [5], p. 114) 

F{x) = -, — cxD < a; < 00, ci > 0. 

ci + cxpl-xj 

(iv) If n = 2fc and a = fe 7^ 0, then Theorem 1 implies 

Xk,2k + <i = Xk+i^2k ~ <2 iff F{x) = -— -— , C2 > 0. 

1 + C2 exp{— a;/KJ 

The next result addresses the case when ([3]) involves non-adjacent (two spac- 
ings away) order statistics. 

Theorem 2 Let k be a fixed integer such that 1 < fc < n — 2 and a > 

and b > be two fixed real numbers such that {a,b) ^ (0,0). Assume X is a 
random variable with continuous cdf F(x) and S,i and ^2 o.i'e independent stan- 
dard exponential variables, which are also independent from Xk^n cind Xk+2,n, 
respectively. Then 

Xk,n + aCl = -'^fc+2,n - bS,2 (6) 



holds if and only if the quintile function Q{x) satisfies 

Q{x)^\og\cx^^{l-x)-'^''-''-^A~W{x;k,n), < x < 1, (7) 

where c > is an arbitrary constant and 

{ d\og{{n - 2k + l)x + k + 1} /oi._Li 

n ZK + i ^g^ 

{a + b)kx, if n^2k + l, 

where d = bk{n ~ k) + a{n — fc — l)(fc + 1). 

Remark. Note that if n = 2k + 1 and a = b — 1/fc, then Theorem 2 impUes 

Xk,2k+i + T" = Xk+2.2k+i - T- iff Qi^) = log \ "i r ~ 2a;. 

k k |_ 1 — a; J 

A number of characterizations of the logistic distribution use distributional 
relations between X and order statistics with positive and negative exponential 
random shifts. George and Mudholkar [4] (see also [8]) proved that X is standard 
logistic if and only if 

X = Xi,2 + e or X = ^2,2 - e- 

More generally, Lin and Hu 8 established that, under some smoothness condi- 
tions, X is standard logistic if and only if 

n — 1 ^ 

X = Xi n + / ^, 
1 

3 = 1 ■' 

where S,j for j = 1,2, . . . ,n — 1 are independent standard exponential variables, 
which are also independent from Xin- Ahsanullah et al. [1] proved that the 
standard logistic distribution is characterized by the following distributional 
equality holding for a fixed k, such that 1 < fc < n — 1, 



x^x,„ + ^iiA^-^ 



.=1 ^ pi ^ 

where Xi^n, ^i(j)j a-nd ^2(j) for j — 1,2, . . . ,n — 1 are mutually independent 
and all ^'s are standard exponential. Finally, Zykov and Nevzorov ^lOj obtained 
characterizations based on cither 

^ = ^«,n-C or X + ^^X„,n- (9) 

Exploring the so-called _F"-scheme (cf. [6], Lecture 25), we will study char- 
acterizations of logistic and related distributions based on one extension of ([9|) . 



Theorem 3 Let a > and b > be two fixed real numbers such that 
(a, &) 7^ (0)0) and F(x) be a continuous distribution function. Suppose Yi 
and Yi are independent random variables with distribution functions F°'{x) and 
F^{x), respectively for a > and /3 > 0. Furthermore, let ^i and ^2 be inde- 
pendent standard exponential variables, which are also independent from Yi and 
m.ax{Yi,Y2}, respectively. Then 

Fi + <i = maxlFi , ^2 } - b^2 (10) 

holds if and only if the quintile function Q{x) satisfies 

Q{x)^\og{cx^°'{l-xf^)-'^], 0<2:<1, (11) 

where c > is an arbitrary constant and d = [a{a + P) + ba]/ /3. 
Remarks, (i) Setting a = I and b — in (fTO)) . yields the cdf 

/ r /3 n^/^ 

F(x) = 1 — exp < -(x — c) > : c < x < 00. 

V I " + /S // 

(ii) li ba = P and d = 1, then (|lip is the quintile function of the Type I 
generalized logistic cdf (cf. 5 , p. 140) 

r (x) ~ 7 , — cxD < a; < 00, C2 > 0. 

\C2 +exp{-a;}y 

The special case of Theorem 3 when a = 1 and /3 = n — 1 is particularly 
noteworthy. 

Corollary Let a > and b > be two real numbers such that (a, b) ^ (0, 0). 
Assume X is a random variable with continuous cdf F{x) and ^1 and ^2 are 
independent standard exponential variables, which are also independent from X 
and Xn,n, respectively. Then 

X + a^i = Xn,n - b^2, n>2, (12) 

holds if and only if the quintile function Q{x) satisfies 

g(x)=log{cx''(l-x"-i)-'*}, 0<a;<l, (13) 

where c > is an arbitrary constant and d = {b + na)/{n — 1). 

The rest of the paper is organized as follows. In Section 2, we present an 
auxiliary result before proving Theorem 1. Section 3 deals with the case of 
non-adjacent order statistics. The proof of Theorem 3 is given in Section 4. We 
summarize the findings and discuss potential future work in the last section. 



2 Adjacent order statistics 

Let Z\ and Z2 be two random variables with continuous distribution functions 
G and iJ, respectively. Let ^1 and ^2 be independent standard exponential 
variables, which are also independent from Z\ and Z2, respectively. 

Lemma Let a and b be two real numbers such that (a, b) ^ (0, 0). Then 

Zi + a^i = Z2 - b^2, (14) 

if and only if 

bG'{x) + aH'{x) = G{x) - H{x). (15) 



Proof. Let us first prove that (iMt implies ^T5\j . Assume both a 7^ and 
6 7^ 0. It is not difficult to see that ((T4)) is equivalent to 

— / G{u) exp < > du ^ — / H{u) exp < — - — > du, 

a J -00 I a J bj^ I J 

which, in turn, implies 

/ G(u) exp < — > du = a exp < > / JT[u) exp < > du. 

J -00 ^a) [ ab j J^ lb) 

Since G and H are continuous, we can differentiate with respect to x both sides 
of the above equation and obtain 

exp {-|} [bG{x) + aH(x)] = ^^ H H{u) exp {-^} du. 

Differentiating with respect to x again (using the continuity of H{x)) we obtain 
(jlSp . The proof when exactly one of a and b equals zero is similar and is omitted. 
To prove that (ITS]) implies dH]), one needs to follow the steps in the above 
proof in reverse order. The proof of the lemma is complete. 

Further on we denote Fi,n{x) = P{Xi^n < x) for 1 < i < n. Recall that 

^»,n(^)=E('')^'(^)(l-^(^))""'' l<^<n. (16) 

Proof of Theorem 1. First, we will prove that Q implies ([5]). Assume 
both a ^ {) and b ^ Q. It follows from Lemma 1 with G{x) = Fk^n{x) and 
H{x) = Fk+i^n{x), making use of ([TC|. that 

&F;„(a:)+a^^+i,„(a;) = Fk^^) ~ Fk+i^x) (17) 

'fjF''ix)il^Fix)r~\ 



On the other hand, for the left-hand side of ([T7)) . using again P^ . we obtain 
bF'{x)+aF;,+,,^{x) (18) 



(a + b) 



dx 



j=k+i 



^ r)F^{x){i-F{x)r-^ 



dx 



F''{x){l~F{x)y 



i — k 



(a + 6)nf" ^)f''{x){1~F{x))''-''-'F'{x) 



n — 1 
fc-1 



(1-n^)) 



n — 1 



^(xl 



+6nF'=-i(x)(l - F{x))''-''^^F'{x) 
an\''',^^F''{x){l-F{x))''-''-'F'{x)+bn(''^2^)F''-\x){l-F{x))"'^F'{x) 



F''-'{x){l - F{x))''-''-'F'{x)[a{n - k)F{x) + bk{l ~ F{x))]. 



Equating the right-hand sides of (fT7|) and (fT8|) . and dividing throughout by 
(^)F'=(x)(l - F{x)Y-^ ^ 0, we obtain 

bk—^ + a(n - k) ^-- = 1, 

F(a;) ^ ^ 1 - i^(2:) 

which, upon integration with respect to x, yields 
log {F^^a;)(l - i^lx))-"'"-*^) } 



X — c. 



(19) 



for an arbitrary constant c > 0. Finally, replacing a; with Q{x) in the last 
equation, we obtain ([5|). The proof for the cases when either a or 6 is zero is 
similar and is omitted here. 

To prove the "if part of the theorem, i.e., ([S]) implies (|3]), one needs to start 
with (|T9l) and repeat the steps of the above proof going backwards to (flT)) . 



3 Non-adjacent order statistics 

In the beginning we will prove that ([6]) yields ([7]). Assume both a 7^ and 
6 7^ 0. Referring to (fT6| . we have 



Fk,n{x) - Fk+2,n{x) 



(20) 



1]f\x){1 - F{x)r-^ +(kll )^'^'(^)(1 - F{x)r-^'^ 



On the other hand, similarly to (fT8)) . we obtain 



bF'{x) + aFl+^,^{x) 



bk[ ', ]F^-\x){l - F{x)y^F'{x) 



(21) 



a{n-k-l)[ '" ]F^+\x){l^F{x)Y-^-'^F'{x). 



Lemma 1 with G{x) — Fk.n{x) and H{x) = Fk+2.n{x) implies that the left- 
hand sides of (|20|) and (|2T|) are equal. Equating the right-hand sides of these 
equations, it is not difficult to see that 

^ (fc + l)F(x)(l - Fix)f + {n- k)F\x){l - F{x)) 
^ ' bk{k + 1)(1 - F{x)y + a{n-k- l){n - k)F^{x) ' 

Replacing x with Q{x) and using that F'{Q{x)) = 1/Q'{x), we obtain 

_ a{n - fc - l){n - k)x^ + bkjk + 1)(1 - xf 
^^'~ x{l-x){k + l + {n-2k-l)x) ■ ^ ' 

If n ^ 2k + 1, then ((22|) can be written as 

,, . a(n — fc — 1) bk a{n — k — l){k + 1) + bk{n — k) 
^ ^""^ " 1^ + ~ fc + 1 + (n - 2fc - l)a: ' ^ ^ 

After integrating with respect to x we obtain ([7]) with the upper branch of (j8]). 
In the case n = 2k + 1, (|22|) becomes 



Q'(a;) ^ -^ + ^ _ (a + 6)fc (24) 

and Q follows after integrating with respect to x again. The proof when exactly 
one of a and b equals zero is omitted. 

Starting with (|23l) and (|24|) and following the steps of the above proof going 
backwards, one can verify that ([7]) yields (jH]). 



4 F" scheme 

Let us first prove that (fTO|) implies pT|) . Assume a ^ and & ^ 0. Furthermore, 
assume /3 = 1. It follows from Lemma 1 with G{x) — F"-{x) and H{x) = 
F°'+^{x) that 

baF°"^{x)F'{x) + a{a + l)F°'{x)F' {x) = F"(a;) - F"+i(x), 

which, upon division by F"^^{x) ^ 0, becomes 

F'{x) [ba + a{a + l)F(a;)] = F{x){l - F{x)). 

This, for F{x){l — F{x)) ^ 0, can be written as 

1 ba (a -|- 6)a + a 

F'(x) ~ F{x) ^ 1 - F{x) 

Replacing x with the quintile function Q{x), we obtain 

Q' OcH — + ^^^ ■ 

X 1 — X 



H 



Integrating both sides with respect to x, we have 

Q{x) = log [cx'"'{l - a;)-(-+b)"--| , (25) 

where c > is an arbitrary constant. This completes the proof when /3 = 
1. To extend this proof to the general case /3 > 0, consider the independent 
random variables ¥{ and F2 with cdf's (iJ(x))"/^ and H{x), respectively, where 
H{x^) = F{x). It follows from (|25|) that the quintile function H~^{x) satisfies 

'i(x) = log {cx''"/^(l - a;)-('^+'')"//3-'^| . (26) 

Now, taking into account that Q{x) = H~^{x^), we see that (P5|) yields ([TT|) . 
The proof of the cases when either a and b is zero is omitted. 

It is not difficult to verify that (TTTI) implies ^TO\i , following the steps of the 
above proof in reverse order. 

5 Concluding remarks 

We obtained characterizations based on distributional equalities between two 
order statistics plus or minus multiples of independent standard exponential 
variables. The resulting family of distributions includes as its members expo- 
nential, logistic, and generalized logistic. In case of adjacent order statistics, 
the quintile function of the underlying distribution takes on a compact explicit 
form. If the order statistics are two spacings away, the expression for the quin- 
tile function includes an additional additive term, which is linear if the order 
statistics are on both sides of the median. Further calculations show that if 
higher-order spacings are involved then the resulting quintile function will have 
more additive terms. 

Exploring the so called i^"-scheme, we studied the distribution of a contin- 
uous X for which X plus a multiple of a standard exponential variable equals 
the maximum order statistic minus a multiple of another standard exponential 
variable. The obtained results generalize those of Zykov and Nevzorov (2010). 

One area of future work will be to study in more detail the distributions 
with quintile functions (O, Q, and (ITT]) . It also remains to be seen what 
the distribution of X is if the corresponding order statistics differ by a linear 
combination of standard exponential variables. Some results in this direction 
are given in Ahsanullah et al. T. 
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